Mon. Not. R. Astron. Soc. 000, 000-000 (2001) Printed 1 February 2008 (MN I*TeX style file vl.4) 



Cosmological constraints from the masses and abundances 
of galaxies 



U. Seljak 



Department of Physics, Princeton University, Princeton, NJ 



, USA 



(N 
O 
O 
(N 



> 



o 

O 

:^ 

Oh 

6 



X 



1 February 2008 



ABSTRACT 

We place limits on the mean density of the universe and the effective slope of 
the linear power spectrum around a megaparsec scale nefi by comparing the universal 
mass function to the observed luminosity function. Numerical simulations suggest that 
the dark matter halo mass function at small scales depends only on ^m{ncs + 3) inde- 
pendent of the overall power spectrum normalization. Matching the halo abundance 
to the observed luminosity function requires knowing the relation between the virial 
mass and luminosity (separately for early and late type galaxies) and the fraction of 
galaxies that reside in larger halos such as groups and clusters, all of which can be 
extracted from the galaxy-galaxy lensing. We apply the recently derived values from 
SDSS and find ri„j(nofr -I- 3) = (0.15 ± 0.05)/(l — /dh), where /dh accounts for the pos- 
sibility that some fraction of halos may be dark or without bright central galaxy. A 
model with Vl„i = 0.25 and primordial n — 0.8 or with fi™ — 0.2 and n — 1 agrees well 
with these constraints even in the absence of dark halos, although with the current 
data somewhat higher values for flm and n are also acceptable. 



1 INTRODUCTION 



Halo mass function has been long recognized as a power- 



ful probe of cosmology since the seminal work by Press & 



Schechter (1974). Most of the applications so far have fo- 



cused on clusters, which are easy to detect in X-rays and for 
which the observed X-ray temperatures correlate well with 
the cluster mass. While there is currently some uncertainty 
in the normalization of this relation, upcoming X-ray and 
weak lensing observations should provide empirical means 
to calibrate it. Since the clusters are the most massive halos 
formed in the universe they lie on the exponential tail of 
the mass function, whose amplitude depends mainly on the 
overall normalization of the power spectrum and the den- 
sity parameter Qm- On the other hand, for the halo mass 
below the nonlinear mass the mass function depends only 
on the density of the universe and on the effective slope of 
the linear power spectrum at that sc ale and is independent of 



lites inside a larger halo and one must correct for that, since 
the halo mass function only counts isolated halos. Estimat- 
ing this fraction is difficult for small halos. 

In the mass ranee 10" - IO^^A/q, typical for L* galax- 
ies, the luminosity function is well determined and theo- 
retical models suggest that each of these halos should host 
a bright galaxy at the center (e.g. Kauffmann et al. 199t , 
Benson et al. 2000). The main challenge in this range is to 



the power spectrum normalization ( Press fc Schechter 1974 



Sheth fc Tormen 1999, Jenkins et al. 2001) 



To constrain these cosmological parameters with mass 
functioij i one must be able to determine masses and abun - 



determine the relation between a galaxy luminosity and a 
halo mass and the fraction of these galaxies that belong to 
larger halos such as groups and clusters. Even though we 
have good dynamical probes of mass within the optical re- 
gion of a galaxy, such as TuUy-Fisher relation for late types 
and Faber-Jackson relation or strong lensing for early types, 
the relation between optical masses and virial masses is more 
uncertain, since the virial mass depends on the adopted den- 
sity profile and a typical virial radius is a factor of 10 larger 
than optical radius. Virial masses thus cannot be observa- 
tionally constrained from the optical or HI studies alone, 
which was used in previous work on the mass function deter- 
mina tion on galactic scales (Sonzalez et al. 2000 



dances hf halns below the cluster scale rirniips with masses 

between 10^^ — 10^*Mq are difficult to observe directly, since 
they contain only a handful of galaxies, and their abundance 
is quite uncertain in this range. Below 10^^ Mq many of the 
halos may not host a bright galaxy because of effects that 
prevent either gas cooling or star formation, such as UV 
background radiation, feedback or insufficient surface den- 
sity for star formation. Moreover, a significant fraction of 
galaxies corresponding to these halo masses may be satel- 



Kochanek 



2001| ). 

An alternative approach adopted here is to use virial 
masses derived directly from the galaxy-galaxy (g-g) lensing. 
In this method one uses tangential distortions of background 
galaxies by the foreground galaxy to place limits on the 
mass distributions around galaxies. A recent SDSS analysis 
presents galaxy-galaxy lensing results using a much larger 
sample than previously available, allowing a detailed study 
of the relation between mass and light for several luminosity 
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bands and morphological types ( McKay et al. 2001 ). The 
latter is particularly important since the relation between 
mass and luminosity differs significantly between early and 
late type galaxies and one must extract the relations sepa- 
rately before combining them together. The data are most 
sensitive to 100-200/i~'^kpc transverse separations, which is 



a typic al virial radius of a 10^^ — W^^h ^ Mf^ halo. Second biased relative t o fainter ones (e.g. Norberg et al. 2001 



mean by co mparing the relation between galaxy light and 
convergence (Wilson et al. 2001). However, if galaxies are a 
biased tracer of dark matter then only a combination flm/b 
can be determined and if 6 > 1 this will underestimate flm- 
Since this is a luminosity weighted statistic it will give a 
larger weight to brighter galaxies, which are known to be 



complication is the fraction of galaxies that are not in iso- 



lated halos. This can also be extracted from g-g lensing: if 
there is a significant lensing signal around the galaxies at 
separations above 200ft~^kpc then it signals a presence o f 
groups and clusters around them (Guzik & Seljak 2002). 



This allows one to determine the fraction of galaxies of a 
given luminosity in these larger halos. Together thus g-g 
lensing provides all the necessary information needed for a 
quantitative study of halo mass function on galactic scales. 

It is worth comparing this approach to the traditional 
mass to light rati o (M/L) method to obtain the density pa - 
rameter fim (e.g. Bahcall et al. 1995, [Carlberg et al. 1997 ). 
This often assumes that M/L extracted from some type of 
objects, such as groups or clusters, can be applied to the 
global M/L. Since we can measure the total luminosity den- 
sity in the universe we can obtain the total matter density 
by multiplying the two. However, M/L depends on both the 
halo mass and scale on which one measures it. While light is 
concentrated to the inner parts of the clusters mass contin- 
ues to i ncrease, so M/L for anv given object increases with 



radius. [Even adopting the virial masses, defined so that only 
baryons within that radius can condense to make stars, there 
is no reason why Mvir/i should not depend on L. Theoreti- 
cal models predict M^ir/L to be at the minimum at galactic 
mass es, where cooling and s tar formation are m ost efficient 



(e.g. Kauffmann et al. 199!: , Benson et al. 200C ) . It is possi 



ble t hat there are many small halos that have no light at all 
(e.i 



Jimenez et al. 1997), which is why we cannot see them, 
but which contribute to the mass of the universe. Moreover, 
there may be mass in the universe that is not associated 
with any collapsed structures at all. For current generation 
of simulations only about 40% of the total mass has bee n 
resolved into halos above 10^^/i~^Mq (lenkins et al. 2001). 
It is not obvious that as the resolution of simulations in- 
creases this fraction converges to unity, which is what is 
commo nly assumed in the forms of ma ss function (Press & 



Schechter 1974 



Sheth 



Tormen 199£| ), since some fraction 

of the mass could remain in a diffuse form. All this makes 
any extrapolation of A/vir/i to the total luminosity density 
as a method to deduce the density of the universe highly 
uncertain. One can attempt to correct for this using simula- 
tions (e.g. Bahcall et al. 2000[), but this relies on the ability 



of these to reproduce the light distribution in the universe 
across a large dynamic range of masses and scales. This is 
one venue to pursue in the future as simulations and mod- 
elling of physical processes improve. 

The alternative way is to measure mass and light over 
a much larger volume, so that it becomes representative for 
the whole universe. Currently there are no reliable meth- 
ods that can measure the mean mass directly on such large 
scales. The closest method to achieve this goal is gravita- 
tional lensing, which measures ellipticity distortions of back- 
ground galaxies. This method however cannot measure the 
mean component of the matter, which does not produce 
shear. So instead one must look at fluctuations around the 



havi et al. 2001) 



Since it is difficult to detect a signal on large scales with 
gravitational lensing, where linear biasing applies, one must 
compare the fluctuations in light and mass at smaller scales. 
At any given scale halos of a c ertain mass dominate the fluc- 
tuations (see e.g. Seljak 2000 for a discussion of this in the 



context of halo models) . For example, on a megaparsec scale 
the dominant fluctuations come from groups and clusters, 
while on 50kpc scale the galactic halos dominate the fluc- 
tuations. This means that even if the galaxies are unbiased 
the relation between light and mass will be appropriate only 
for the halos that dominate the fluctuations at that smooth- 
ing scale. As discussed above this mass to light relation may 
not be the universal one, since most of the mass may be in 
smaller objects or in diffuse structures, which have too weak 
signal to be detected through the weak lensing fluctuations. 
There is no preferred scale at which one should evaluate 
M/L to multiply it with the total galaxy light to derive the 
mean density of the universe. If M/L increa ses with halo 



mass above L,, as suggeste d by observations ( Gtirardi et al 
200l|, iGuzik fc Seljak 2002| ), then the convergence-light cor- 



relation function will be more extended than that of light 
itself and there is some observational evidence for this effect 



(Wilson et al. 2001). It is not clear however whether this 



leads to an underestimate or overestimate of Q,m, since the 
trend of M/L to increase with mass is likely to be reversed 
below L* and a lot mass associated with either diffuse struc- 
tures or small halos may not be associated with any light at 
all. 

Above arguments suggest that global mass to light ra- 
tio method of determining the density of the universe cannot 
be derived without making additional assumptions. In the 
approach presented here we instead limit the analysis only 
to the galaxies around L, which are well studied. We com- 
bine the virial mass to luminosity relations extracted from 
the SDSS data with the SDSS luminosity function of L* 
galaxies and compare that to the universal mass function 
to place limits on cosmological models. The analysis is done 
entirely within the SDSS data set, which reduces the un- 
certainties related to the photometric calibrations and color 
transformations, which usually plague luminosity function 
comparisons. 



2 RELATION BETWEEN GALAXY AND 
HALO ABUNDANCES 

The halo mass function describes the number density of ha- 
los as a function of mass. It can be written as 



dluM M''^ ' dlnM 



(1) 



where p is the mean matter density of the universe, M is 
the virial mass of the halo and n{M) is the spatial number 
density of halos of a given mass M. We introduced a func- 
tion /(cr), which has a universal form independent of the 
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Figure 1. n^ff +3 for models with = 0.3, f2ij = 0.04, h = 0.7, 
n = 1 (solid, top), n = 0.9 (solid, middle) and = 0.25, 

= 0.04, h = 0.65, n = 0.8 (solid, bottom). All the trans- 
fer functions were computed using CMBFAST. Also shown is 
f[a{M)] (dashed), which is essentially 0.3 over this range of 



power spectrum, matter density, normalization or redshift if 
written as a function of rms variance of linear density field, 



a^{M)=4n / P {k)W R{k)k'^ dk . 



(2) 



Here Wnik) is the Fourier transform of the spherical top hat 
window with radius R, chosen such that it encloses the mass 
M — AtvR^p/Z and P(k) is the linear power spectrum. 

The relation between mass and rms variance depends 
on the Unear power spectrum. For a smooth power spec- 
trum one can locally approximate it as P{k) oc k""^^ and 
the relation is 



dlnM 



^cff + 3 
6 ■ 



(3) 



For CDM models on galactic scales the effective slope rioff 
ranges between -1.5 and -2.5 (figure |l|). 

The universality of the mass fun ction has been recently 
investigated by a n umber of authors ( |Sheth fc Tormen 1999[ 
Jenkins et al. 2001, White 2001). It has been shown that the 



mass function is indeed universal for a broad range of cosmo- 



logical nodels. Jenkins et al. (2001) propose the following 
form, 



/(a) = 0.315 exp[ 



Ina"' -h0.61|^'**l 



(4) 



which they argue is universal if mass is expressed in terms 
of the virial radius where overdensity is 200 in units of mean 
density. It is remarkable that the mass function in this form 
is almost constant for all halos with Incr > 0. Here we are 
interested in halos on galactic scales where Ino" ~ 1, so this 
limit a pplies and we can take / ^ 0.3, a universal value 



which varies only weakly with halo mass, as shown in figure 
|l] (dashed line). 



To relate the theoretical halo abundance with the ob- 
served galaxy abundance we must assume a relation between 
the galaxy luminosity and its halo mass. Here we use a di- 
rect probe of the halo vi rial mass obtained f rom the galaxy- 
galaxy lensing of SDSS ( |McKay et al. 200l| ). Detailed mod- 
elling of CDM profiles shows that these observations are 
best fitted with L, galaxies (which dominate the luminos- 
ity distributi on of that sample) ha ving a virial mass around 



lO^^/i"^M0 ( iGuzik fc Seljak 2002| ). The virial mass strongly 



depends on the morphology or color: for a given luminosity 
early type galaxies can be significantly more massive than 
late type galaxies. The differences are reduced in red bands 
(r', i' and z'), where we focus our analysis as well, but they 
can still vary by a factor of 2-3. This is to a large extent a 
consequence of stellar age, which modifies the stellar mass to 
light ratio by a similar factor, with early type galaxies having 
a higher stellar mass to hght ratio than late type galaxies. 
Current data are consistent with stellar mass to virial mass 
being independent of stellar age (or morphology) . Neverthe- 
less, in the absence of stellar mass information for individual 
galaxies, one must extract all the relations separately for the 
two types before adding the m together. For exam ple, for an 



L* = 2.1 X IO^^^^^Lq in i' (|Blanton et al. 200lD one finds 



virial M2oonm ~ 4.3 x 10^^ /i~ Mq for late type galaxies and 
Af20onm = 1.2 X lQ^^h~^ Mq for early type galaxies, where 
we have corrected for a 25% increase from M200 to M2oon™ 
assuming Q,m = 0.3 fiducial model. 

An alternative and more model dependent estimate of 
virial masses come s from the TuUy-Fishe r (TF) relation for 
late type galaxies ( Gliovanelli et al. 1 9971) or Faber-Jackso n 



(FJ) relation for early type galaxies (Bernardi et al. 2001). 
These observe galaxy properties in the inner lQh~ kpc and 
do not directly measure the virial mass of the halo in which 
the galaxy sits. At these radii the dynamical effect of baryons 
on the rotation curves (as well as on the dark matter dis- 
tribution itself) is important and modifies the relation be- 
tween the rotation velocity and virial mass. One can how- 
ever model this assuming stellar mass to light ratio and dark 
matter profile. The first comes from the stellar population 
synth esis models and is rather uncertain because of age and 
IMF (Bruzual & Chariot 1993), while the latter is obtained 
from the cosmological dark matter sim ulations and depend s 
on the assumed cosmological model (Bullock et al. 2001). 



In addition, the response of dark matter to baryon cooling 
must be included and is often modelled with adiabatic con- 



traction (Blumenthal et al. 1986). Adopting standard values 
for the dark matter profiles and the stellar mass to light ra- 
tios one finds that at L* the rotation velocity decreases by 
1.8 from optical to virial radius both for early and lat e type 
galax ies, in good agreement with g-g lensing results (Seljak 



2002| ) 

Galaxy-galaxy lensing can also be used to determine the 
slope P of the relation between mass and luminosity, 



M* " VIT/ 



(5) 



where M» is the mass associated with L« galaxy. Using SDSS 
data one finds /3 = 1.4 ± 0.2 in red bands, which becomes 
= 1.2 ± 0.2 after correcting for a lumin osity dependent 



fract ion of early type galaxies in the sample (Guzik & Seljak 



2002). This is valid only for early type galaxies between L» 
and 7L,, which dominate the g-g lensing signal. Note that 
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this differs significantly from FJ relation, which would pre- 
dict /3 ~ 2/3, but agrees well with the d etailed mode lling 
of rotation curves for early type galaxies (Seljak 2002). To 
determine (3 around L» for late type galaxies we cannot use 
g-g lensing, since the signal is too weak to be detected as 
a function of lur ninosity. Instead we u se TF relation which 
gives Li qc v^p^ ( Giovanelli et al. 19971 ). We model the con- 
tributions from stellar disk and adiabatic response of dark 
matter to disk formation to relate between the rotation ve- 
locity at the optical radius and virial velocity (or mass) . We 
assume NFW profile with C2oo ~ 10 and use T/ — 1.5h, 
which were shown to reproduc e well the vi rial velocity con- 



straint from g-g lensing at L, (Seljak 2002). We find (3 



around L*, which is the value we adopt below. 

The third parameter that we need is the fraction 7 of 
galaxies that are at the centers of isolated galactic halos, as 
opposed to larger halos such as groups and clusters. Uni- 
versal mass functions from N-body simulations only count 
isolated halos, so one must correct for the fraction of galaxies 
at a given luminosity that are in larger groups and clusters. 
This fraction can be determined from the relative contribu- 
tion to galajcy-galaxy lensing at small and large separations. 
Above 200-300/i~^kpc the signal is dominated by groups and 
clusters, while below it is dominated by individual galac- 
tic halos. Analysis of SDSS galaxy-galaxy lensing data finds 
7 = 0.72 ±0.1 for early type g alaxies and 7 = 0.93 ±0.1 
for late type galaxies around L, (Guzik & Seljak 200S). The 



error includes various systematic uncertanties, of which the 
galaxy occupation as a function of group and cluster mass is 
the most important. Thus about 10-30% of L, galaxies re- 
side in groups and clusters as defined in N-body simulations 
and one must reduce the galaxy abundance by this fraction 
when relating it to the halo abundance. 

The abundance of galaxies of a given luminosity can be 
extracted from the luminosity function dn{L)/dln L, which 
determines the abundance per logarithmic interval of lumi- 
nosity. It is often fitted to the Schechter form. 



dn 
dlnL 



(r) 



a + l 



exp(— L/i,). 



(6) 



For L, galaxies the abundance is "I>(L*) — 4'*/e, where 
e = 2.718 is the natural logarithmic base constant. From 
the early SDSS analysis the values for 0* are 1.46 x 10~^, 
1.28 X 10"^ and 1.27 x 10"^ (with a 10% error) in units of 
/I'^Mpc"'^ for r', i' and z' , respectively, while the val ues for 



a are around -1.2 to -1.25 with a few percent error ( Blan- 



ton et eI. 2001). We choose these 3 bands because they show 



least variation in virial mass to light ratio between early and 
late type galaxies, minimizing the color dependence of the 
signal. In addition, the redshift evolution corrections in red 
bands are smaller than in g' or u' . These have been sug- 
gested as one possible reason for the discrepancy between 
the 2d F lu minosity function in b.j and SDSS in g' (Norberg 
et al. 2lp0l|). 

The last ingredient that is needed is the fraction of early 
type (5°) and late type (£') galaxies a s a function of magni 



tude. As shown in Strateva et al. (2001) this fraction depends 
on luminosity, so that the early type galaxies dominate at 
the bright end and the late type galaxies dominate at the 
faint end. This statement is only valid for red bands and a 
reverse trend is seen in u' , while g' shows comparable frac- 
tions almost independent of luminosity. In red bands around 



L, the fraction is somewhat higher for the early types. While 
the transition between the early and late types is not well 
defined the transitional types (SO-Sa) do not dominate the 
counts, so here we will simply assume a bimodal early/late 
distribution (^"^ ± ^'=1), rather than attempt to model the 
whole range of stellar ages and morphologies. This could be 
improved in the future as larger statistical samples are ob- 
tained and is particularly important for spirals, which have 
a larger scatter in the stellar ages. Early type galaxies are 
more homogeneous and old, so their scatter in stellar mass 
to light ratio should be small. 

We can now combine the above equations to relate halo 
and galaxy abundances. This gives 



dn 



dn 



7<E>(L) _ n,ff ± 3 



(7) 



dlnM (5d\nL (5 6 M2oon„. 

Using p = f7„pcrit=2.77 x IO^/i^MqMpc"^ one finds the 
minimum mean density at a given luminosity is given by 



^m{ncH ± 3) = 0.3 



<I>(L)M2oon„ 



fl3 \ IA X IOIo/iWqMpc- 



(8) 



Note that one must add up the contribution from both early 
and late type galaxies separately, where the two contribu- 
tions have to be evaluated at equal mass, not luminosity. If 
there are dark halos without a bright galaxy at the center 
the above expression become inequality and one can only 
place a lower limit on ^2m(?^off ± 3). We parametrize this 
uncertainty with the fraction of dark halos /dh, which in 
general is a function of halo mass. 

We can evaluate this expression at several different val- 
ues for halo mass. The virial mass of an early type galaxy at 
L, = 2.1 X 10^°fe~^Ln in i' is M2oon„ = 1.2 x lO^^h'^MQ 
( [Guzik fc Seljak 200!^ ). Using = 1.2, 7° = 0.72 and 
= 0.6 one finds ni,{n^a ± 3) = 0.14/(1 - /dh). To this we 
must add the contribution from late types at the same mass. 
At L, their virial mass is M2oon„ = 4.3 X lO^/i-^M©. Us- 
ing = 1 one finds that for M2oon™ = 1-2 x IO^^/i'^Mq 
the corresponding luminosity is 3L* . At this luminosity the 
fraction of late type galaxies in the sample is already small, 
^' ~ 0.2. In addition, from the luminosity function in equa- 
tion (^ one can see the abundance of 3Lt galaxies has de- 
creased by a factor of 10 relative to that of L» . So late type 
galaxies do not actually add much to the limit above and 
together we find f2ni(?icff ± 3) = 0.15/(1 — /dh)- These con- 
straints are evaluated in i' , but one finds similar constraints 
also in r' and z' . The error budget is dominated by the er- 
rors on Af20on„ , 7 and which combined give about 30% 
uncertainty. 

We can repeat the same analysis one magnitude above 
and below L*. At L = 2.5L* = 5.2 x 10^°/i"^Lq the sample 
is dominated by early type galaxies, so ^'^ ~ 0.8, P'^ ~ 1.2 
and M = 3.6 X 10^^/i"^Mq. The fraction of these galaxies 
in isolated galactic halos i s not well determined, but is likely 
to be larger than at L, (Guzik & Seljak 2002), so we will 
assume 7*^ = 0.9. This gives r2in(?ieft ± 3) > 0.13/(1 — /dh), 
where we have ignored the very small contribution from the 
late type galaxies. The effective slope is about 5-10% higher 
than at 1.2 x lO^^/i'^M© (fi gure uX), so the obtained value 
is about 20% lower than the value obtained at L*. If we 
assumed 7" does not differ from that at L* we find the two 
estimates are in perfect agreement. This is quite impressive 
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given that the masses and abundances change by a factor of 
several. The error is comparable to the error at L«. 

One magnitude below L« the sample is dominated by 



late type galaxies, for which we use f3 



1, M = 1.7 X 



10"/i"^Mq and 7' ~ 0.9. Adopting 5' = 0.8 leads to 
f^m(woff + 3) = 0.11. To this we have to add the contribution 
from early types at M = 1.7 x lO"/i"^M0. If /3'= = 1.2 ex- 
tends to this mass range this mass corresponds to L ~ O.IL, 
and at this luminosity the fraction of early type galaxies is 
about 10%. This increases the above estimate by 20%, so 
^m{naS + 3) = 0.13/(1 — fdh)- This estimate is more uncer- 
tain, since both /3''° and 7''° have not been directly measured 
over this range. In addition, the large scatter in stellar ages 
for late type spirals leads to a scatter in mass to luminosity 
relation. Note that at this mass one expects rics -I- 3 to de- 
crease by about 20% relative to 1.2 x IO^^/i'^Mq (figure [i|), 
so this constraint is actually very similar to the one at L, 
based on early type galaxies, even though the masses differ 
by almost an order of magnitude. 

From the above analysis we find that over the range 
of masses between 1.6 x lO" /i'^Mq to 3.6 x 10^^/i~^Mq 
the cosmological constraints on Qm/{1 — f^h) are all very 
similar, 



n„,(neff +3)(l-/dh) =0.15±0.05. 



(9) 



The good agreement found over a wide range of mass sug- 
gests that the shape of the mass function agrees well with 
the one predicted from cosmological simulations, assuming 
the fraction of dark halos /dh, if different from zero, does 
not vary over this mass range. 



3 DISCUSSION 

In this paper we propose a method to derive cosmologi- 
cal constraints from the virial masses and abundance of 
L, galaxies and apply it to early SDSS observations. The 
method differs from other methods using mass to light ratio 
M/L in that it only uses this information around L, galax- 
ies and not the overall luminosity density. This sidesteps the 
uncertainties related to the variation of M/L with luminos- 
ity L. The main observational inputs are relation between 
virial mass and luminosity around L, as a function of mor- 
phological type and the fraction of these galaxies in isolated 
halos as opposed to groups and clusters, all of which can 
be extracted from g-g lensing. Another essential ingredient 
is the luminosity function of galaxies around L* as a func- 
tion of morphological type, which can be obtained from the 
same data as g-g lensing information. The main theoreti- 
cal input is the halo mass function, whic h has been shown 



to be universal by a number of studie s ( Sheth & Tormen 



1999| , [Jenkins et al. 200l| , [White 200 1| ). The abundance of 
halos depends only on the density parameter flm and the 
effective slope of the linear power spectrum n^s through the 
combination Qm{neti -1-3). 

In the absence of dark halos the obtained constraint 
^lm{ncH + 3) = 0.15 ± 0.5 is low compared to the predic- 
tions of ACDM model with flm ~ 0.3, h = 0.7 and n — 1, 
which gives f2,„(ncff -I- 3) = 0.28 at these scales. This is ex- 
cluded by the current constraints, unless a significant frac- 
tion of halos is dark over this range of masses. To bring the 
models into a better agreement with this constraint one can 



either lower the effective slope or the mean density. The for- 
mer can be lowered by reducing the primordial spectrum 
slope n, decreasing the shape parameter F (which depends 
on fim ~ ^ dm + ^b. and Hubble constant Ho, see e.g. 



Eisenstein fc Hu 1990 or introducing warm dark matter 



iodo, Ostriker, 



Turok 2001). For example, we find that 
a model with Qm = 0.25, h = 0.65, Qb = 0.04 and n = 0.8 
gives a good agreement with the observational constraints, 
but somewhat lower values of Qm ~ 0.2 with n — 0.9 — 1.0 
are also acceptable. For warm dark matter models, which 
suppress power on sma ll scales, we compute tran sfer func- 
tions using CMBFAST ( [Seljak fc Zaldarriaga 199^ ) and find 
that the effective slope at the galactic mass scale remains al- 
most unchanged for > 500eV and drops to rics — —2.5 
at — 250eV. Such lo w masses are probably excluded 
from Ly — a forest studies (Narayanan et al. 200C), although 
this must be confirmed with a more careful error analysis of 
Ly ~ a forest constraints. 

There are several possible sources of uncertainty that 
can affect the current constraints derived above. First there 
is the possibility that the virial masses used here are too low. 
This is certainly possible for late type galaxies, which have a 
weak signal in g-g lensing and for which a factor of 2 increase 
in mass is less than a 2-a excursion. For early type galaxies 
the statistical error on the virial mass is 20%, so a factor of 
2 excursion is unlikely, although some systematic uncertain- 
ties remain in the g-g lensing analysis. The agreement be- 
tween early and late type galaxies implies that the mass scale 
must be changed for both types. Furthermore, any increase 
in the virial mass would affect the agreement between g-g 
lensing masses and TuUy-Fisher or Faber- Jackson relation. 
A change in mass by a factor of 2 would reduce the optical 
rotatio n velocity to virial rotation velocity ratio from 1.8 
to 1.4 (Seljak 2002). Within the context of adiabatic com- 
pression models such a ratio is likely to be too small to be 
explained with CDM profiles and stellar mass to light ratios 
as expected from stellar population synthesis models with 
reasonable IMF. This possibility would thus require one to 
give up a successful prediction of CDM models, that of the 
structure of dark matter halos in the outer parts. 

Another possibility is that the relation between the halo 
mass and luminosity M oc L'' is shallower than assumed 
here, (3 ~ 0.6. This would be the case if virial mass to light 
ratio was decreasing with luminosity. Both g-g lensing and 
modelling of optical relations are consistent with /3 = 1 — 1.5 
at L, and above and theoretical models also predict /3 > 1 i n 



this regime (e.g. Kauffmann et al. 1999, Benson et al. 2000) 



While M/L probably decreases with L well below L,, this 
is unlikely to be the case over the range of luminosities of 
interest here. If = 1 instead of 1.2 used here this would 
increase the lower limit on f2m(noff -I- 3) by 20%. It seems 
similarly unlikely that the luminosity function can be wrong 
by more than 30%, unless there is a large fraction of low sur- 
face brightness galaxies that are missed by the SDSS survey. 
There are still important calibration differences between dif- 
ferent surveys, whi ch can cause a misrn atch in derived lu- 
minosity functions ( Norberg et al. 2001 ), but these are not 
relevant for our analysis since we derive the g-g lensing re- 
lation between light and luminosity using the same sample 
that is used for luminosity function as well. 

Another uncertainty is the division into early and late 
type galaxies and their associated fractions, which is some- 
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what artificial, since there is a continuous range of colors and 
light concentrations in the sample. What is really relevant 
is the stellar mass of the galaxy, which seems to correlate 
quite well with the virial mass. For late type galaxies the 
range of stellar ages is large and this introduces a spread in 
the luminosity for a given stellar mass, so representing them 
with a single mass to light ratio may not be accurate. This 
is less of an issue for early type galaxies, which are very old 
and for which the spread in the stellar mass to light ratio is 
small. On the other hand, early type galaxies tend to reside 
in denser enviroments and the fraction of these in groups 
and clusters is larger. This correction is also somewhat un- 
certain, since one must assume how the groups an d clusters 
are po pulat ed with these galaxies to determine it (Guzik & 
Seljak '. 002). There is also the issue whether the mass pro- 



file of galaxies within groups and clusters differs from that 
of the same luminosity in the fie ld within inner lOOkpc. Ob- 
servations (Guzik fc Seljak 2002) and numerical simulations 



probably hosts several fainter galaxies rather than a single 
bright galaxy at the center. It is difficult to determine the 
abundance of such groups directly from the optical data, 
since it is not obvious which groups are virialized to satisfy 
the halo definition in an N-body simulation and which are 
just a collection of galaxies approaching each other for the 
first time (as for example the local group). The fact that 
the constraints derived here agree from 1.6 X lO"/i"^M0 
to 3.6 X 10^^ Mq suggests that over this range of masses 
the simplest possibility with /dh = is consistent with the 
data. 

The prospects to determine the fraction of dark halos 
directly seem difficult. The only direct way to observe these 
is through the gravitational lensing, but in the absence of 
a significant b aryon condensation such halos are inefficient 
strong lenses (Kochanek & White 2001). The halos that do 



( Ghigna et al. 2000] ) suggest they do not differ much, but 
the uncertainties are still large. The extreme case is when 
there is actually no mass attached to individual galaxies in- 
side groups and clusters. Then one should use 7*^ = 1, which 
would increase the estimated Slm(?icff + 3) by 30%. 

The choice of the virial radius and its associated mass, 
defined differently by different authors, has a minor effect 
on our results. From the simulations it is usually defined 
by friends-of-friends algorithm (assuming a specific value of 
linking parameter, e.g. 0.2), but it is not always clear how 
this relates to the specific r nean density within the virial 
radius. For consistency with Jenkins et al. (2001) we define 
the virial radius as the value where the mean overdensity 
is (5 = 200, while the observed masses are usually expressed 
in terms of overdensity relative to the critical density. For a 
given h alo profile one can convert between the two, but the 



cool and form a disk without making stars could be some- 
what more efficient for a given halo mass, but are expected 
to have lower halo masses. Most of the lenses are bright early 
type galaxies which both reside in massive halos and have a 
significant baryonic contribution to the lensing cross-section. 
It is thus possible that even if some fraction of halos around 
10^^ Mq exists without a central galaxy they may not 
have been detected so far with strong lensing. Current sur- 
veys such as SDSS may provide better limits on the fraction 
of dark halos as a function of halo mass. 

It is clear from the above discussion that the errors in 
the current analysis are still quite large, although the fact 
that the constraints are consistent over a range of masses 
increases the confidence level of the final result. It is inter- 
esting that the constraints obtained are in a good agreement 
with the cluster gas fraction determina tion of the matter 
density (Erdogdu, Ettori, & Lahav 2002) and with t he red- 



convers .on depends on the assumed density parameter 
Fortun ately the differences are not very important on galac- 



shift distor tions and bias determination from 2dF (Verde 



et al. 2001). They are also compar able or somewhat h i gher 



than those usin g global M/L ratio ( Bahcall et al. 200C , Wil 



tic scal3s, where halos are highly concentrated and where 
the mass around the virial radius only slowly grows with 
radius (logarithmically in the limit of a large concentration 



son et al. 2001), although this method, as discussed above. 



may well be biased. Similarly, a tilted CDM model m ay help 



aleviate some of the small scale problems with CDM (Alam 



where 1he slope at virial radius approaches -3). For exam- 



Bullock, & Weinberg 2001). This perhaps indicates the need 



pie, the difference in mass between S = 20057m for 0,m. ~ 0.4 
and Q,n = 0.2 in a halo of c = 12 is 10%, so this is not the 
dominating source of error. This effect has been included in 
the estimates above. 

The remaining uncertainty is the fraction of dark halos 
in the universe. Theoretically one would not expect halos to 
be dark around i*, where the efficiency of cooling and star 
formation is high. This is supported by the fact that for the 
halos that we do see a large fraction of the baryons within 



the viri al radius has converted into stars (Guzik & Seljak 
2002|)~Kt lower halo masses, below 10"/i 'Af0, the frac- 
tion of dark halos may increase because some of the formed 
disks r nay be stable against star formation ( Verde, Oh , & 
Jimene: 



2002). At higher halo masses above 



10"^" 



'Mr, 



© 



one enters the group regime, cooling becomes less efficient, 
and a significant fraction of these halos may not host a bright 
galaxy at the center. For example, applying the same anal- 
ysis as in this paper to TL* early type galaxies one finds the 
abundance of halos to be several times below that expected 
from the mass function, a consequence of exponentially de- 
creasing luminosity function at the bright end. This is not 
surprising, since a large fraction of the groups in this range 



for a somewhat lower Q,m or rioft than previously suggested 
Qm = 0.3 - 0.4, n = 1 model. 

While the error from the method presented here is still 
large, the prospects to improve it are good. Currently the 
errors are dominated by observationally determined param- 
eters M2oon„, 4>*, P and 7. These errors are dominated by 
statistics and were obtained by using only 5% of the final 
SDSS sample. With the full sample one can reduce the error 
on each of these significantly, as well as extend the observ- 
able range to the lower luminosity galaxies. With the full 
sample one can also study the morphological dependence 
of the mass-luminosity relation in more detail, which as we 
have shown here plays an important role in the analysis. 
With spectroscopic information it should be able to extract 
stellar mass information for each galaxy separately and use 
g-g lensing to relate the stellar mass to the virial mass di- 
rectly without splitting the sample into morphological types 
as done here. Theoretical errors are mainly caused by the 
accuracy of the mass function over this mass range, but the 
uncertainty is already at a 10% level and can be further 
improved with simulations that cover a wider range of cos- 
mological models. The remaining uncertainty is the fraction 
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of dark halos, which should also be determined with a com- 
parable accuracy (or shown to be negligible). In this case 
the method presented here may become an accurate test 
of matter density and slope of the power spectrum on one 
mcgaparsec scale. 
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